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The pattern transition induced by lattice anisotropy (LA) and magnetic impurities is computa- 
tionally observed in near-optimally doped d-wave superconductors (DSCs). For the single impurity 
case, a transition from the checkerboard to stripe pattern can be induced even with a very weak 
LA. Moreover, the modulation period of eight lattice constants (8a) in the spin order coincides with 



'TIS ■ neutron scattering data. For the two-impurity case, an orientation transition from the longitudinal 

O . impurity-pinned stripe into the transverse pattern is observed when the LA ratio reaches some crit- 

ical value. At the critical point, it is found that the structures around magnetic impurities could 
^ ■ restore checkerboard patterns. These results indicate that the formation of stripes in DSCs might 

cn : 

CN ■ induced by various effects, and could be tunable experimentally. 
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The inhomogeneous phases in unconventional superconductors have attracted much at- 



tention recently. Various experiments reportec 



the presence of stripe or checkerboard mod- 



12|. Neutron scattering (NS) measurements 



ulations in copper oxide-based compounds [l 
on cuprates such as La2_x'-?/Ndj^Sr^.Cu04 (LNSCO), La2_xSrxCu04 (LSCO), La2_xBax.CuO 
(LBCO), and Y-Ba-Cu-0 (YBCO) family indicated the existence of incommensurate mag- 
netic peaks, which has led to discussions of the existence of a stripe phase jl -3|. Several scan- 
ning tunneling microscopy (STM) experiments observed checkerboard-like charge-density 
wave (CDW) modulations with a period of roiighly four lattice constants (4a). Experimen- 
tally, both dispersive [5] and non-dispersive Q-qj modulated patterns have been observed in 
cuprates. The corresponding patterns were proposed to be understood in terms of quasipar- 



ticle interference (QPI) 13|, or in terms of static or fluctuating stripes 15|, respectively. 



Currently, the issue of the nature of these modulated patterns is still under debate. 

It is known that lattice anisotropy (LA) are ubiquitous in various cuprates. For in- 
stance, strong a-b axis asymmetry of both the normal and superconducting state electronic 



properties have been observed in the YBCO family 



I6I-II8]. Up to now, the relationship 



among lattice distortions, incommensurability and stripes has been intensively studied in 



cuprates 



19|-|26|. Based on the anisotropic Hubbard (t^ 7^ ty) and — ty — — Jy models. 



Normand and Kampf et al. considered LA as a possible origin of the stripe formation in 



the cuprate superconductors 



19| . Moreover, Becca et al. have shown that stripes and spin 



incommensurabilities are favored by LA [22]. In this work, we explore the effect of LA 
on the patterns around magnetic impurities in d-wave superconductors (DSCs). Including 
the competition and coexistence among the DSC, spin, and charge orders, we study the 
system by self-consistently solving the Bogoliubov - de Gennes (BdG) equations based on 
an anisotropic Hubbard-type model. 

We start from the two dimensional tx—ty — t' — U — V model, which consists of two parts, 
H = Ho + Himp- The Hamiltonian Hq and Himp describe the superconductor and magnetic 
impurities, respectively, which is given by 

Ho = ~^tij{clcj^ + H.c.) 

ijcr 



Himp = ^heff{i){cl^Ci^~Cyil), 



(1) 



Here Qcr annihilates an electron of spin a at the ith site. The hopping integral tij takes 
tx or ty between nearest neighbor (NN) pairs along x or y direction, and t' between next- 
nearest neighbor (NNN) pairs. U is the on-site Coulomb repulsion interaction, /i is the 
chemical potential, which is determined self-consistently in the calculation. The local ef- 
fective field heff{i) is introduced to model the exchange coupling between conducting elec- 
trons and the impurity spin, where we have treated the impurity spin as a Ising-like one. 
Some similar model had been employed to study the effects of magnetic impurities on 
cuprate superconductors 



states well 



271 . |28|, which can qualitatively explain the observed impurity 
29|. Therefore, we employed the above model in this work. Experimentally, the 
ratio of the lattice constants along x and y directions is 6/a ~ 1.01 (Ref. [30]), thus the 
corresponding effect.ve hopping .ntegrals and are also anisotropic. Aeeord.ng to the 

local-density approximation band calculation |31l |. the ratio of hopping integrals is estimated 
as tx/ty ~ {h/aY ~ 1.04. In this work, we consider effects of LA on the patterns around mag- 
netic impurities by tuning the ratio rj = ty/tx, which is near the above estimated value, as 
shown below. The self-consistent mean- field parameters are given by rii = < cJ^Qa- >, 
the magnetization = (1/2) (< c|^Cj^ > — < cj^Cii DSC order parameter 

Aij = {V/2) < Ci^Cji — Ci^Cj^ > with V the phenomenological pairing interaction. 
The Hamiltonian H can be diagonalized by solving the following BdG equations. 



= i?^. (2) 



where the quasiparticle wave function \E'j = (wj^, ^4)^. The spin-dependent single-particle 
Hamihonian reads Hi^^ = -txSi+xj - ^^i+y,j - t'^i+r'j + [Y.i„, (^heff{i)6i^i^ + Urii^^ - ^\5iy 
Here the subscripts x and y denote the unit vector directing along x- or y-direction NN 
bonds, t' denotes the unit vector directing along four NNN bonds, and im is the position of 
the impurity site. The self-consistent parameters are given by rii-^ = |u^p/(i?„), = 
E„ - fiE^)], and A,,- = ^ EnK^^t + v^{u]^]tanhi^), where f{E) = 1/(1 + e^^) 

is the Fermi-Dirac distribution function. Hereafter, the length is measured in units of the 
lattice constant a, and the energy in units oit = tx- The pairing interaction is chosen a.sV = 
1.0 to guarantee that the superconducting order Aq — 0.08t, comparable with the observed 
Tc in cuprate superconductors. The on-site Coulomb repulsion U takes value 2.44, which is 
the effective value near optimal doping in the hubbard-type model 13|, |32j. We study LA 



effect on the patterns around impurities in DSCs near optimal doping with the filling factor 



Uf = J2ia'^ia'^io-/{NxNy) = 0.85 (i.e., the hole doping x = 0.15), where A^^, A^^ are the hnear 
dimension of the unit cell. The BdG equations are solved self-consistently for a square lattice 
of 24 X 24 sites, and the periodic boundary conditions are adopted. The numerical calculation 
is performed at a very low temperature, T = lO^^K, to extract the low-energy physics. The 
local effective field is taken to be h^ff at the impurity site and zero otherwise. The DSC 
order parameter at the zth site is defined as = (Ai^i+e.^ + ^^i^i-e^ — ^i,i+ey — Aj^j_ej^)/4, 
and the spin order parameter is Mj = (— l)*mj. 

In FIG.[I1 we plot the spatial distributions of the DSC, spin, and charge orders around the 
impurity site. As one can see [FIG.[T](a), (c), and (e)], all the three orders exhibit checker- 
board modulations around the magnetic impurity. Similar to the nonmagnetic impurity case, 
a SDW with a period of 8a checkerboard pattern is observed around the magnetic impurity, 
which is in agreement with the NS experiment j^. However, one notes that the checkerboard 
pattern of the DSC order can also be induced by the magnetic impurity, and a weak associ- 
ated CDW pattern is observed, which is different from the nonmagnetic impurity case |l3 |. 
Moreover, the modulated DSC and CDW orders share the same periodicity 4a. On the 
whole, the DSC order has a local minimum at the impurity site while the amplitudes of the 
CDW and SDW orders reach global maxima. This is the common feature of orders around 
the magnetic impurity despite various parameters. Therefore, in view of these features, we 
clearly see the relationship of competition and coexistence between antiferromagnetic and 
DSC orderings. 




Figure 1: (color online) The surface plots of orders around the magnetic impurity on a unit cell of 
size 24 X 24 sites with and without lattice anisotropy. (a), (c), and (e) are the spatial distributions 
of the DSC, spin and charge orders without LA. (b), (d), and (f) are the same plots but the lattice 
anisotropy ratio ry = ty/tx = 0.99. Here we take U = 2.44, x = 0.15, t' = —0.25, and /ig// = 3. 
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Below we discuss the effect of LA on the pattern around a single impurity [See Fig. [T] (b), 
(d), and (f)]. As one can see, a symmetry-broken transition from the checkerboard to stripe 
pattern is induced. In this case, even if there is a very weak LA, this transition takes place. 
The reason lies in the fact that LA ty < t,j. breaks the symmetry of the system, which leads to 
the formation of stripes along y direction favorable. Therefore, we see that LA is favorable to 
the formation of stripes in DSCs. In addition, one notes that the anisotropy-induced SDW 
stripes also show a modulation with periodicity 8a, which coexists with the DSC order. This 
indicates that the 8a modulations appearing in the SDW order is a robust feature in spite 
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of different parameters. This conclusion is consistent with our previous study 

Now let us turn to consider the two-impurity case. Without LA, an x-orientated stripe 
structure can be induced due to the pinning effect of the magnetic impurities which have 
been placed along x direction [Not shown but similar to FIG. [2] (a), (c), and (e)]. Simi- 
lar to the single impurity case, the impurity-pinned SDW stripe also show a modulation 
with periodicity 8a, while the coexisting DSC and CDW stripes share the same periodicity 
4a. For the case of two magnetic impurities with LA present, one could expect that the 
impurity-pinned effect would compete with LA effect, leading to orientation transition from 
the longitudinal (x direction) impurity-pinned stripe into the transverse (y direction) pat- 
tern when the LA ratio reaches some critical value. This speculation is confirmed by the 
following calculations, and the critical value is found to be rj* ~ 0.982, which is close to the 
above estimated value ty/tx ~ {h/a)^^ ~ 0.9615. 

Accordingly, the ratio of lattice constants is estimated as 6/a ~ 1.005. Thus, we see that 
only a weak LA is needed to observe this pattern transition in the disorder-pinned stripe 




Figure 2: (color onhne) The surface plots of orders around two magnetic impurities in the presence 
of LA. (a), (c), and (e) are plots in the case with -q = 0.983, while (b), (d), and (f) with r] = 0.981. 
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y. If adopting the estimated value by Normand and Kampf 



33] . one gets r] = ty/t^ = \cos{tt — 2<l>)| ~ 0.9848, with an 



angular distortion angle $ = 5°. As shown in FIG. [2] (a), (c) and (e), when the LA ratio is 
larger than the critical value (a weak anisotropy), t] = 0.983 > r]*, the x-orientated stripe 
structure is still stable. This case is very similar to the case without LA, except that the 
modulation amplitudes are slightly larger than those in the latter case. However, for the 
relatively strong anisotropy with rj = 0.981 < rj*, the orientation transition from x-stripe to 
y-stripe pattern takes place, as shown in FIG. [2] (b), (d), and (f). Clearly, in this case, LA 
effect dominates over the impurity pinning effect, and the y-stripe pattern is energetically 
more favorable. Also, as shown by all those stripe plots, the checkerboard-like modulations 
are still visible, especially in the DSC and SDW orders, which originate from the effect of 
QPI. Thus, we see that one can produce different modulated patterns by inserting magnetic 
impurities or tuning the LA in DSCs while the modulation periodicity holds the same. 

For comparison, we also consider the cases with various strength of effective field, heff = 1 
or 10, and with a long distance between the two impurities. In these cases, we obtain similar 
results, and reach the same conclusion as discussed above, except that the critical value t]* 
changes more or less. For the cases with h(,ff = 1 and with a 8a distance between the two 
impurities, we happen to get the exact critical values of the LA ratio, i.e., rj* = 0.984 and 
0.998, which are both larger than the previous case with heff = 3. Meanwhile, for the case 
with heff = 10, the corresponding critical values is found to be t]* ~ 0.981, which is smaller 
than the above cases. The reason lies in the fact that for a weaker effective field or a longer 
distance, the impurity pinning effect is reduced, while an opposite tendency is produced by 
a stronger effective field. 

At the exact critical value of the LA, it is found that the structures around magnetic 
impurities restore checkerboard patterns because the competition between the impurity pin- 
ning effect and LA reaches the balance point, as shown in FIG. [31 In this case, it seems 
that the two impurities could induce their own checkerboard patterns independently, with- 
out interference with each other. These results indicate that one could observe the stripe- 
checkerboard-stripe transition under some conditions by tuning rj in impurity-substituted 
DSCs. 

In summery, we have studied the effect of LA on the patterns around magnetic impuri- 
ties based on the tx — ty — t' — U — V model. It is demonstrated that LA could induce pattern 
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Figure 3: (color online) The surface plots of orders around two magnetic impurities at the critical 
LA ratio with different distances between the two impurities, (a), (c), and (e) are plots in the case 
with r]* = 0.984, while (b), (d), and (f) the case with r/* = 0.998. 

transition around impurities in near-optimally doped DSCs. In the single impurity case, it 
is found that even a very weak LA could induce a transition from the checkerboard to stripe 
pattern, because the symmetry of the system is broken. Modulated SDW pattern with 8a 
periodicity is observed, which coincides with the NS data. Meanwhile, the modulated DSC 
and CDW orders share the same periodicity 4a. For the two-impurity case, a transition 
from the x-direction impurity-pinned stripe into the y-direction stripe pattern is observed 
as the LA ratio reaches the critical value rj*. At the exact critical values, it is found that the 
structures around magnetic impurities could restore checkerboard patterns. It is expected 
that these phenomena could be observed in the STM and NS experiments on DSCs under 
varying pressure or temperature. 

This work is supported by the National Nature Science Foundation of China (Grants No. 
10904063 and No. 10804047). CDC would also like to thank the 973 Project (Project No. 
2006CB601002). 
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